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Abstract 

We are interested in the homogenization of elastic-electric coupling 
equation, with rapidly oscillating coefficients, in periodically perforated 
piezoelectric body. We justify the two first terms in the usual asymptotic 
development of the problem solution. For the main convergence results 
of this paper, we use the notion of two-scale convergence. A two-scale 
homogenized system is obtained as the limit of the periodic problem. 
While in the static limit the method provides homogenized electroelastic 
coefficients whicht coincide with those deduced from other homogenization 
techniques (asymptotic homogenization 0, F-convergence [16j). 
Key words. Homogenization; Piezoelectricity; Perforations 



1 Introduction 

Composites and perforated (lattice) materials are widely used in many practi- 
cal applications, such as aircraft, civil engineering, electrotechnics, and many 
others. These materials are with a large number of heterogeneities (inclusions 
or holes), and in strong contrast to continum materials, their behavior is defini- 
tively influenced by micromechanical events. 

The first goal of this work we study the homogenization of the equation of 
the elastic-electric coupling with rapidly oscillating coefficients in a periodically 
perforated domain. The homogenized of this problem for a fixed domain has al- 
ready been studied, by the author (Feng and Wu.^Oj, Castillero and aZZ.[5j, Ruan 
and all. [14j). But in this work we give new convergence results concerning the 
same model by using homogenization technique of " two-scale convergence ", 
which permits us to conclude the limit problem, the approximation of final state 
is altrered by a constant named as the volum fraction which depends on the pro- 
portion of material in the perforated domain and is equal to f when there are 
no holes. 
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The second goal of this paper is to apply the technique of a formal asymptotic 
homogenization, to determine the effective elastic, piezoelectric and dielectric 
modulous of periodic medium. The final formulae for the effective parameters 
are given in a relatively simpler closed form. 

The third goal of this paper is to also establish a corrector-type theorem, 

which permits to replace the sequence by its "two-scale" limit using the result 
on the strong convergence, and permits to justify the two first terms in the usual 
asymptotic expansion of the solution. In the last section we treat of the enregy 
aspect of our problem. 

2 Homogenization problem 

Throughout this paper in the Sobolev space of real- valued fimctions that 

are measurable and square summable in with respect of the Lebesgue measure. 
We denote by (Y) the space of infinitely differentiable functions in that 
are periodic of Y. Then, L'^{Y) (respectively, H^{Y)) is the completion for the 
norm of ^^(y) (respectively, H\Y)) of Cj°°(y). 

2.1 Geometric of the medium 

Let n c be a bounded three dimensional domain with the boundary F = dSl. 

We denote x the macroscopic variable and by y = |- the microscopic variable. 
Let us define fie of periodically perforated subdomains of a bounded open set 
O. The period of 0^ is eY*, where Y* is a subset of the unit cube Y = (0, 1)^, 
which represents the solid or material domain, S* obtaindcd by K-periodicity 
from Y*, is a smooth connected (the material is in one piece) open set in R^. 
Denoting by xiv) the caracteristic function of S* ( Y-periodic) , in well be 
defined analj^tically by 



2.2 Model problem 

We adopt the convention of Einstein for the summation of repeated indices, we 
use Latin indices, understood from 1 to 3, we note by the fields of displace- 
ment in elastic, and by ip'^ of electric potentiel. The equations of equilibrium 
and Gauss's law of electrostatics in the absence of free charges, written as 



f2e = {xeO, x(-) = i} 




(1) 




(2) 
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where f G L (rig). The second-order stress tensor cr^ = i^ij)' elec- 
tric displacement vector D"^ = (-Df), are linearly related to the second-order 
strain tensor Ski{u) = ^{di-ui + diUk) and the electric field vector dkf^ by the 
constitutive law 

<ji-^',f') = ctjkiM^') + e%,jdkf' in n^, 
i<i,j,k,l<3, 

where (div a^y = d.,af^, div D'^ = d,Dl 9, = .t = (a;,) € fl. And the 
material proprieties are given by the fourth-order stiffness tensor c^^f.^ measured 
at constant electric field, the elastic coefScients satisfy the following symmetries 
and ellipticity uniformily in e, 



^ijkli^) ~ Cijkl{x, g), 

^ijkl ^ '^'jikl ~ ^klij ~ '^ijlk^ (a\ 

3ac a.c.[e) > : cfj,.iXijXki > acXijXij, VXij = Xji e M. 

The third-order piezoelectric tensor ef^^, (the coupled tensor), verify the follow- 
ing symmetry, 

^ijk ~ ^ijkix, - ), 

eijk{x,y)eL°°{Cl;Ct{Y)). 

The second-order electric tensor dfj (dielectric permittivity), measured at con- 
stant strain, verify the conditions of symmetric and ellipticity uniformily by 
£, 



(6) 



di,(x,y)eL-(fl;Cj(F)), 

Baa 7^ ad{s) > : dfjXiXj > aaXiXi, VXi e K. 

2.3 Variational problem 

Introducing the two Hilbert spaces 

Ve(Oe) = { V e Hi(f2e), V = on d^l} 

We(Oe) = { ^ G H^{n^), ^ = on dn} 

With two norms : || . ||v^(n^) = || . ||Hi(n), || • ||we(a,) = || • The 
variational problem is defined by : 

Find (u^<^^) G Vei^s) x We{Cls), such that 

(7) 

ae((u^ <^^), (v, ^)) = Le(v, V) V(v, ^) G Ve(n.) X W,{Q,), 
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where 



It is pointed out that under assumptions (jH)-©-®, the variational problem 
(O-® have a unique solution (u^,y)^) £ Ve(rie) x We(rie), corresponding the 
saddle point of this functional (see [12j ) : 



(v,^)-^2/ (c'(v,v) + 2e^(u,V)-#(^,V)) da:- / f v dx, 



where 



c^(u,v) = cfjfei 5fe;(v) 
e''(u,'(/') = efj.; Sfcz(u) SjV' 
d'^{i;,'ip) ^ dlj dii} djil> 

2.4 A priori estimates 

In order to prove the main convergence results of this paper we use the notion of 
two-scale convergence which was introduced in [10] and developed further in [1] . 
The idea of this convergence is based in first step by taking a priori estimates 
for displacement field and the electric potentiel. The second step we use the 
relatively compact property with the classical procedure of prolongation (wich 
is the extension by from fi^ to fi). Finally we pass the limit e — > 0, in order 
to obtain the homogenized and the local problems in same time. 



Proposition 1 By using the two equivalent norms of\'^{^l^), We(f2e), for any 
sequence of solution (u'^,(/?^)e C V£(f2e) x Ws;{^e) of variational problem 
(0). Then, this solution is hounded, and we have this a priori estimate uniformly 
by e 

II I|hi(o.) + II \\HHn,)< C, (9) 
where C is constant strictement positive and independent by e. 

Proof. 

By choosing v = u'^ and ip = Lp^ in variational formulae ([T])-®, and by using 
the Korn's and Poincare's inequalities in perforated domains (see Oleinik et al. 
[13] for the Korn's inequality and Allaire- Murat [2l for Poincare's inequality), 
we see that and Lp^ are bounded, by a constant which does not depend on e. 
For other details see [12j . 
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3 Two-scale convergence 



We denote by ~ the extension by zero in the holes fl ~ fl^. The sequence of 
solution (u'',(/3^)e C Ve(17e) X Wei^le) of variational problem ©-([l]) verify ^ 
and, in this case, by adding the relatively compact property and elementary 
properties of two-scale convergence, imply 

Lemma 1 1. There exists u{x) £ Hj(rj) and (p{x) S Hq{^) such that, the 

two sequences (u^)e, two-scale converge to x(2/)u(a;), x(y)iy9(a;), re- 

spectively. 

2. There exists Ui{x,y) G L2[rj; Hj(r*)/K], ipi{x,y) e L^[n; HI{Y*)/R] 
such that, 

V u*^ — > x(2/)[VxU(x) + Vj,Ui(a;, y)] in two-scale sense 
\/ ip^ —>■ x{y)[^x^ix) + S7y(pi{x, y)] in two-scale sense 

3. We have 

s (u^) x(2/)[s^(u(a;)) + Sy{ui{x,y))] in two-scale sense 
index x or y means that the derivatives are with respect to the variable. 
Proof. For details see [J, [TD], [U and [H]. 

Corollary 1 The sequence (u'^)e>o {resp. {lp^)^-^^), converge weakly to a limit 
6u {resp. 6ip) in L^(r2) {resp. L^(r2)). 

Remark 1 Let p e L'^{Y), define p^{x) = pif), and (w'^)e C L'^{Q) two-scale 
converge to a limit v G L^{fl x Y). Then {p^v^)^ two-scale converges to a limit 
pv (see U^). 

From last results, we can state next theorem 

~ ^ ~ ~ 

Theorem 1 The sequences (u'^)e, {s (u'^))^, (</5'^)e and (V '^')e two-scale con- 
verge to x{y)u{x), x{y)[sx{u) + Sy{ui)], x{y)v>{x) and xivWx'P + "^y^fii] re- 
spectively, where (u(x), Ui(a;,y), ipi{x,y)) are the unique solutions in 
B.l{n) X L2[f7;Hj(y*)/M] x H^{Q) x L^p; Hl{Y*)/R] of the following two- 
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scale homogenized system; 
d f 



_d_ 

dxi Jy* 



{-eiki{x,y)[ski,xiu) + Ski,y{ui)] + dij{x, y)[dj.x<p + dj^yipi\}dy] = m f2, 



d 

--^{cijki{x,y)[skLx{vL) + Sfc;,y(ui)] + ekij{x,y)[dk,xV + dk,yVi]} = mVLxY*, 



d 

-K—{-e-iki{x,y)[ski,x{^) + Ski,y{ni)] + dij{x,y)[dj^x^ + d^^yipi]} = in x F*, 

(10) 

and we have this boundary conditions 

\i{x) = on dil, 

ip{x) — on dfl, 

{ctjki{x,y)[ski,x{u) + Ski,y{ui)]+ekij{x,y)[dk,xV + 9k.yVi]}-'nj = on dY* - dY, 

{-eiki{x,y)[ski,x{vi) + Ski,y{\ii)]+dij{x,y)[dj^x'P + dj,yipi]}.ni = on dY* - dY. 

(11) 

y Ui(a;, y) is Y — periodic, 

(12) 

y ^ ipi{x,y) is Y — periodic, 

where 6 is the volum fraction of material {i.e. 9 =< x >= Jy xiu) dy =\ Y* \), 
I Y I denote mesure ofY. 

The equations (llOp - (fTT1) - P^ are referred to as the two-scale homogenized 
system. 

Proof. From the idea of G.Nguetseng [TU], the test functions in Q-® are 
chossed on the form 

v^(a:) = v(x, -) = v°(a;) +ev^(x, -), 
e £ 

r{x)=i^{x,-)= i;'^{x)+e^\x,-), 
e e 
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where vO G Cff (f7), £ Clf{n), 6 Cff (f7; Cg»(y)) and G Co°°(f]; Cj°°(y)), 
we obtain 

{[Cy7ciSfci(u') +eLj5fc¥'1[s^j,x(v°)(x) + {s,y,j,(vi) +es„-,^(vi)}(x, -)] 



f,{x) [v°{x) + evl{x, -)]dx 

Under the precedent hypotheses, and passing to the two-scale hmit, yields 

[c,jki{x, y)x{y){skLx{u) + Ski,y{ui)) + ektj{x,y)x{y){dk^x'P + dk,y(pi)]x{y) 
[ski,x{v°) + Ski,y{v^)] dx dy 
&iki{x, y)x{y)iski,x{u) + Ski,y{ui)) + d^j{x, y)x{y){dj^x(p + dj,yifii)]x{y) 
[di^xi^^ + di^yip'^] dx dy 

h{x)x{y)v^{x) dx dy. (13) 



n JY 



nJY" 



By definition of x, we have 

[cijkiix,y){ski,x{u) + SkLy{ui)) + ekij{x,y)[dk,xV + 9k,yVi)] 
[sfei,j:(v°) + Sfci,y(v^)] dx dy 
-etki{x, y){ski,x{u) + Ski,y{ui)) + dij{x,y){dj,x(p + dj^yipi)] 
[di,x'>P° + dt^yip'^] dx dy 

Ji{x)vi{x) dx, (14) 



O JY* 



where = Jy x{y) dy, by density of spaces from which we chose the test func- 
tions, the equation (|14p holds true for any v*' G Hj(il), t/j^ G HQ{n), and for 
anyv^ G L2[r2; Hj(r*)/K], V'^ e L^[n; Hl(Y*)/R], 

Integrating by parts, shows that (|14[) is variational formulation associated to 
the two-scale homogenized system 

-dj[ / {cijki{x,y)[ski,x{u) + Ski,y {ui)] + ekij{x,y)[dk,x'f + dk,y(pi} dy] ^ 0fi{x) 



di[ / {~eiki{x,y)[skLx{u) + Ski,y{ui)] + dij{x,y)[dj^x'f + dj.yifi]} dy] =0 



(15) 

We complete (fT5|) by the boundary conditions (fTT|l - (fT2l) . To prove existence 
and uniqueness in (jl4p , by application of the Lax-Milgram lemma, let focus on 
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the coercivity in Hj(f7) x L2[r2; Hj(y*)/R] x H^{n) x L^lfl] Hl{Y*)/M.] of the 
bihnear form defined by the left-hand side of (For a complete demonstration 
see [H]). 



Remark 2 /< is evident that the two-scale homogenized problem if J6<)) - f77)] - IJii)j 
is a system of four equation, four unknown (u, Ui, ip, ipi), each dependent 
on both space variables x and y (i.e. the macroscopic and microscopic scales) 
which are mixed. Although seems to be complicated, it is well-posed system of 
equations. Also it is clear that the two-scale homogenized problem has the same 
form as the original equation. 



The object of new paragraph is to give another form of theorem which is more 
suitable for further physical interpretations. Indeed, we shall eliminate the 
microscopic variable y (one doesn't want to solve the small scale structure), and 
decouple the two-scale homogenized problem (fT0)) - (ITl1) - (|12p in homogenized and 
cell equations. However, it is preferable, from a physical or numerical point of 
view (see [12]). 



4 Derivation of the homogenized coefficients 

Due to the linearity of the original problem, and assuming the regularity in 
variation of the coefficients, we take 



OXfi 

Vi{x,y) = s„,HA^{x))^^\y) + ^^V"(2/), 



(16) 



(17) 



where w™'', (y?", q™'' and -0" are F*-periodic functions in y, independent of x, 
solutions of these two locals problems in Y* 



dyj L 



mh . 



dy 



kl , I mh\\ I I \ 

-Tmh + Sfc;,y(w )J -f efc»j (a;, y) | 



d ( r \ \ kl I rah\\ i / \ 1 

I -e,fci(a;,?/)[T„h + Sfci,j,(w )J +djj(x,?/)-^— I 



= in Y* 



= in Y* 



Y* — periodics. 



(18) 



where 



ih = ^[SkmSih + 5khSim] 1 < k,m,l,h <3. 
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_d_ 



d r 

^1 - eiM{x,y)ski,y{cC) + dij{x,y) 



Skn H 
Sjn 



dip'' 



dyk 
dip' 



9yj 



= in Y* 



in Y* 



Y* — periodics. 



(19) 



However, in general a relation (jl6p - (fT7|) like this does not exist, if we have't 
linearity of problem. 

Now substitue the expansions and ([T7]) in this equation 



--K—{cijkiix,y)[ski,x{u) + SkLy{ui)] + ekij{x,y)[dk,xip + dk,yipi]} = in fixF* 



we obtain 



\c^jki{x,y)\T^h + SkLy{^'^^)\ +efc.u(x,y)-^| 



+ ^^|cijfei(a;,y)sfci,a(q ) + eki]{x,y)\bkn + ^ — J/'" 



(20) 



Calling Tmh the basic of symmetric second order tensors r^^^ — \[5km5ih + 

], where Sij is the Kronecker symbol. 
Analogoulosy, we substitue the expansions (|16p and (|17p in this equation 



d 



dyi 
we obtain 



{-etki {x, y) [skLx (u) + Ski,y (ui)] + (x, y) [c^j^"^ + dj.yipi]} = in x y * , 



' '-| -e,fcz(a;,y)[r^';, + Sfc,,j,(w™'')J +dy(a:,y)^— I 



+ - eiki{x,y)skLyisi ) + dij{x,y)\djn + | 



0. 



(21) 



From the relation ([20|) -(l2i p . after lenghty calculations we arrive at the homog- 
enized (effective) coefficients : 



cfjmh = (c^jki{x, y) [r,^,'^ + sfci,y(w™'')] + efcy(x, y) 



dip 



■mil 



dyk 



Cijki{x,y)ski,y{cC) + ektj{x,y) 



^kn 



dip'' 



dyk 



fSnh = \eiki{x,y)\T^h + Ski,y{v^"'^)\ - d,j{x,y)^—^ 



dfn = ( - e^kl{x,y)skl,y{cC) +dij{x,y) 



dijj" 
9yj 



(22) 
(23) 
(24) 
(25) 
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where (h) = h{y) dy, the measurements on Y* of function h. 

Now we give the results concerning some properties of elasticity homogenized 

tensor. 

Proposition 2 The coefficients of elasticity homogenized tensor = i'^ijki) 
defined by i22\) satisfy: 

a) cfju = c§,j = = 4ik ' VI < J, /c, ; < 3, 

b) There exists aff > 0, such that for all ^, symmetric tensor = 

Proof. 

The part of the symmetry of these coefficients is evident 

ijmh jimh jihm' 

We are interesing in the proof of 

^ijmh mhij 

Following the ideas, we transform the above expression to obtain a symmetric 
form. 

We defined the tensor of second order S, by S*^' = \{ykSi + yiSk), we define 
its defined by 



3x3 iJ'"'' matrix by i?'"'' = Sy(S'^'), it is evident the coefficients of this matrix 



[H'''']jnh = T^h = \ [5km,5lh + SkhSlm] I < k , I , m, k < 3 . 

If we use this new notation, we can rewrite the problem (|18p , as the form given 
as under 

-^{cvki{x,y)ski,y[^""' +^''"') +ekrjix,y)^^} =Oinr*, 

-^{-e^M(x,y)sfei.,(E'"''+w"'') +d,,(x,2/)^} = in F*, ^^^^ 

w'"'', ip"'^ Y* -periodics. 

We introduce the problem functions (w*^ , q'-* ) , solutions of the problem 

d ( ( ; ■ \ d(f>^^ 1 

-^|cfcjQ/3(a;, y)sa/3,y(^i;'-' + w^-'J + eQfei(a;,y)-^^ j = in y*, 

-^{-efea0(x,y)s„;3,,(s^^ +w^^) +4,(x,2/)^} = in F*, ^^^^ 
w*-', (/s'^ F* — periodics. 
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The coefficient of elasticity tensor can be rewritten as 

4mh - J^^ c^jki (x, y)skLy (S"'' + w'"'^) dy + J^^ e^, (x, y) ^ dy. (28) 

Tlie second integral of the right-hand side of precedent expression, is evaluated 
as follows 

ekij{x,y)-^ dy = / ekaf3ix,y)-^ daiSpj dy, 



Y- dyk Jy- ' dVk 

1 f dip"''' ( \ 

2 y ekai3{x,y)—^— \5ai6(3] + SajSpij dy, 

dip"'''' 

ekai3{x,y)— SQ/3,y(w'-') dy 

Y' clyk 

+ / ekM^,y)-^Sai3,y{^'' +^'') dy- (29) 
Jy' oyk 

We use the variationnal formulation of the first equation of problem (pS)) . and 
taking the test function v = w'-' , we obtain 

dp"''' f 

ekai3{x,y)—^ SQ/3,y(w*-') dy ^ Ca0kiix,y)ski,yiT."'''+w"''')safj,y{w^^) dy. 

Y* clyk Jy- 

(30) 

Multiplying the second equation by p"''', and integrating by parts, we have 

f dp"''' f dp"''' dip^-^ 

/ ekai3{x,y)— SQ/3.y(S*-' + w'^) = / dka{x,y)— — dy. (31) 

Jy' oyk Jy* oyk dya 

Regrouping these results, and using the definition ([28]). we derive 



rah 



cfjmh = / c,^kiix,y)skLy{T,"''' +W"''') dy+ I ekrj{x,y)^^ dy. 



Y' JY' 

c,M^,y)ski,y{^""'+^""') dy 



dyk 



Y' 



r dc/j"''' dio^-^ 

Cai3ki{x,y)ski.y{^"''' + w"''')sa,i3,y{vf'^) dy+ / dka{x,y)-^ — dy, 

Jy* oyk dya 



Ca/3fei(x,2;)sfe,.^(I]"" + w'"'')sa^,„(S'^) dy+ / dkc.{x,y)^^^ dy 
Y' Jy* oyk Oya 

Ca,3fei(a;,2;)sfei.y(S'"'' + w™'')s„^,y(w^J) dy, 

Y* 

CaM^,y)skLy(J:""' + w'"'')sap,y{j:'^ +yv'^) dy 



Y* 



Y- 



dp"'''dp'^ 

dka[x,y)—^ ^ — dy. (32) 

dyk dya 
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It is immediate from above, that the coefRcients of elasticity tensor satisfies 

ijmh mhij 

This is the end of the proof of the first section of proposition i.c of symmetry. 
We now study the elhpticity of the coefficients of the elasticity tensor, we recall 
'^ijki is elliptic, if for all the second order tensor Xij symetric (X^ = Xji), we 
have 

3 af > 0, c",,i Xij Xki > a" Xij Xij 
Consider expression ([^ of the tensor cfj^i^, we have 



%mh ^ij ^mh = J Cij„ih{x,y) Xij Xmh dy + J Cijki{x,y)ski,y{^™''') Xij Xmh dy 

Q^rnh 

ektj{x,y)^ Xij Xmh dy, 

oyk 

(x, y) Xij Xmh dy + y)skUy{^'^^Xmh)X,j dy 



ekij[x,y) r Xij dy, 

oyk 

Cijki{x,y)(^SkLy{w) + Pki^Pij dy + J ekij{x,y)-^Pij dy, 

where w = w'^'^Xmh, C = 'p'^'^X^h and Pij = r^mh^mh = Xij . Therfore the 
couple (w, () is a saddle point of the functional J defined by 

(v,V)^ J(v,^) 

dtp dip 
dyi dy^ 



•^(vW^) ^\J^^ {%fci(sy,a(v)+i^y) (sfei,j,(v)+Pfc,)+2eHj^(s.y,y(^^ 
By definition of the saddle point, we have 

J{w,ip) < J{w,() < J{'v,() V(v,-0) periodic functions. 

Or 

J(w,0) = i J Cijki{x,y)(^Sij^y{w) + Pij^ (^Ski,y{w) + Pki^ dy. 
But, if we use the first equation of system (fT8|) . we obtain 

Cijmh Xtj Xmh = 2J(w, C), 

> ^(w,0), 

= ^ y" Cijki{x,y)(^Sij^y{w) + Pij^(^Ski,y{w) + Pki^ dy, 

> 0, 

We have the homogenized elastcity tensor — {cfj^i), which is elliptic. 

Now we give a results concerning some properties of dielectric homogenized 

tensor. 



'0 
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Proposition 3 The coefficients of dielectric homogenized tensor T>^ = [d^] 
defined by \25\) satisfy: 



a) dg=d^,, Vl<z,n<3, 

b) There exists > 0, such that for all vector ^, 

Proof. 

By analgy, we transform these coefficients to obtain a symmetric form 

d" = d" 

The problem p^ . can be rewritten as 

d ( d(yn + i/i") ~| 
--Q;^\^Cijki{x,y)ski,yiq") + ektjix,y) ^— 1 = in Y*, 

-g;^]^- eiki[x,y}skLy(q ) +d,j[x,y) — 1 =Omy, 

q", i/i" Y* - periodics. 

Introducing (q%-0*), in the solution of this local problem 

--Q;i^[ctjki{x,y)sM,y{q') + ektj{x,y) 1 = in F*, 

-^{ - eMx,yWiA^') + = in Y*, ^^^^ 

q*, -0' y* — periodics. 

We can rewrite these coefficients of electric tensor in form given as 

-etkiix,y)ski,y{q"') dy + dij{x,y) ^ dy. (35) 

Jy Jy* oyj 

The first term of the second integral of precedent expression, is evaluated as 
follows 



-eiki{x,y)ski,y{q^) dy = j ~eaki{x,y)ski,y{(f')5ai dy, 
= I -eakiix,y)ski,yiq'')^^ dy 

dya 



eaki{x,y)ski,y{cC)^— dy 



'eaki\x,y)ski,y\q ) ^ dy. 

dya 



(36) 
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Using the variationnal formulation of the second equation of system ([55]) . and 
chossing a test function (p — tp'^ , we obtain 

-eaki{x,y)skLy{q")^^ dy ^ [ daj {x, y) '^^^"^ ^ ^ dy. 



dya Jy' ' dyj dy^ 

Let us now consider the second integral in p6p . Multiplying the first equation 
of system (j34p by and integrating by parts, we have 



-eafc;(x,y)sfc;,^(q") ^^^' "'"^ ^ dy ^ - [ Ckia0ix,y)sa/3,yiq'')ski,yi(p'^) dy. 

Oya JY* 

Finality, we regroup these lasts results, and using the definition (I35p . we obtain 
dfn = I -eiki{x,y)ski,y{q") dy + [ (a:, y) ^ ^ dy, 



dy. 



daj[x,y) ^ — dy - / Ckiai3(x,y)Sai3,y{q jSkLyiq ) dy 



I / J / N "\yn ^ y 1 
+ / d,j{x,y) — dy, 



dyj dya 



da]{x,y) — dy - I Ckiap[x,y)Sai3,y{q )Ski,y[ci ) dy 



dyj dya 
,d{ yn + V"") % 
dyj dya 



+ / daj{x,y) — T^— dy. 



daj[X,y) ^ ^ dy- / CklafS{x,yjSa!3,y(Cl )Skl,y{<i 

oyj clya 



It is clear from that the coefficients of electric tensor is symmetric. 

Now we are interested in the ellipticity of this tensor, recall df^ is elliptic, if for 

all vector Xi, we have 

3 > 0, df^ X, X„ > X, X,. 

We consider the expression (pS)) of tensor d^^, we derive 

df„ X^ Xn ^ J din{x,y) Xi X„ dy - J eiki{x,y)skLyiq") X, X„ dy 

dip" 

dij{x,y)^ — Xi Xn dy, 

dyj 

din{x,y) Xi Xn dy - J eikiix,y)skLy{q"Xn) X, dy 

+ / di.j{x,y) X^dy, 

' dyj 



J d,j{x,y)(Qj + -^"jQi dy ~ J 



eiki{x,y)ski,y{s)Qi dy, 



14 



where ^ = ip^Xn, <^ — f'^Xn and Qi = SinXn = Xi. Or the couple (^,<;) is a 
saddle point of the functionnal G defined by : 



G(v,^) 



(v,V^) ^G(v,7/;) 

2 J {^CtjkiSij^y{v)ski,y{^r)+2ekijSij^y{v)(^Qk+-^^-dij(^Qi+ (^J"*"^) 
By definition of the saddle point, we have 

G(C, ip) < G(C, <r) < G(v, (,) V(v, V) periodic functions. 

Or 

But, if we use the second equation of system ()19p . we have 



dij{x,y){Qj + ^^Qi. ^y^\f^ eiki{x,v)ski,y{s) X^ dy, 

-df„ X.t Xn, 



2 7^. '^-'-V--^ dy, 
1 

< 0. 



We have the dielectric homoginized tensor which is elliptic. 

Now we give the results concerning some properties of piezoelectric homogenized 

tensor. 

Proposition 4 The coefficients o] 'piezoelectric homogenized tensor — (e^^j) 
defined by \23\) satisfy: 



Moreover, we have the identity 



e" = e" 



e" ■ = f" 



Proof. 

By definition of coefficients e^^^ (using the fact that Cijki{x,y) = Cjiki{x,y), 
ekji{x,y) = ekij{x,y)), we have 



e" 



\^c,jki{x,y)skLy{cC') + ekijix,y)(^6kn + } '^V' 

\^Cj^kiix,y)ski,yiq") + ekjiix,y)(^Skn + } '^V^ 



(38) 
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We can rewrite the coefficients , as given as under 



\enij{x,y) ^ ekij{x,y)-^+ CijM{x,y)skUy{c('')^ dy. (39) 



Same, we can rewrite the coefficients f^ij, as form givens under 



\^enij{x,y) + enki{x,y)ski,y{w'^^) - dnt{x,y)-^j dy. (40) 

Using the two variationnals formulations corresponding of problems (jlSp and 
(|19p . and chossing the appropriete test functions, we can directly prove as e^^^ = 

■ 

Finallity, using the three last propositions, we can purpose the altarnative form 
of the principal convergence theorem 
Theorem-Bis {the altarnative form) 

Set (u, (^) solution of the two-scale homogenized problem ilO\) - ni\ )- US\) . then 
(u, if) is defined by that the solution of this homogenized problem 



—div a-^{u,ip) = Of inVt, 

-divD^(u,(^) = inVL, 
where the boundary conditions 

Vl{x) = Q on dft, 

(p{x) — on dfl, 
afj and Df are defined by the homogenized constitutive law 

H I \ H ( \ H 

d 



(41) 



(42) 



(43) 



the homogenized coefficients cf^^i, e!^^j and d^ are defined respectively by (2, 
and mi). 



5 Correctors result 

The corrector results are obtained easily by the two-scale convergence method. 
The objective of the next theorem justify rigorously the two first terms in the 
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usual asymptotic expansion of the solution. 

Following the idea of Allaire [1] , we introduce the following definition 

Definition 1 We call 'ijj{x, y) an admissible test function, if it is Y -periodic, 
and satisfies the following relation 

lim / ^{x.-'fdx^ / / ip{x,y)'^dx dy (44) 



s 



n Jy 



Here we recall the Allaire's lemma 
Lemma 2 (Allaire fljl) : 

Let the function •4>{x,y) G LF'{Q,]C^(Y)), theml:{x,y) is an admissible test func- 
tion in the sense of Definition [7J 

Using this lemma, we obtain the following proposition. 

Proposition 5 The two functions Sij_y{ui(x,y)) and di^yLpi{x,y) are admissi- 
ble test functions in the sense of Definition]^ 

Proof. 

By definition we have 

ui(x,y) = s™;,,,(u(a;))w™'^(2;) + %^q"(y), 

Viix.y) = s^h,.{vL{x))^^'\y) + 



we obtain 



S^j^y{Ul{x,y)) = S„^h,x{u{x))s^j^yiw'^''{y)) + Sij ,y{(f (y)) , 

di^yipi{x,y) = SmhA'^{x))d.,^yip'^^'{y) + ' di^y^p"iy). 

UXn 

Using Lemma[2l Sy.j^(ui(x, y)) and di^yipi{x, y) are the admissible test functions 
in the sense of Definition [TJ 

Theorem 2 We have these two strong convergence results 

{x) — x( — )[u(a;) + Ui(a;, — )] strongly in Hj(r2) 

{x) - x{-)[v{x) + (pi[x,-)\ ^ Q strongly in Hl{n) 
e e 
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Proof. We consider the variational formulation under the following form 



{[cljkiSkiiu')+el,jdk^']s^j{v)^-elkiSki{u')+dljdjip']d,^p}dx = j fi{x)v,{x) dx 

(45) 



Chosing V = u"^ and ip = Lp^ vcl (|45|) . we obtain 

{[4jfe;Sfci(u^)+eLjafc^^]sy(u^)+[-effe;Sfe,(u^)+df^aj^^]a^^^}dT = j f,{x)ul{x) dx 
By simplification, we get 



Mjkii^)st3iu'')ix)ski{u''){x)+d'ij{x)d^ip'' {x)djip'' (x)} dx ^ f,{x)ul{x) dx 

(46) 

By applying we can write 

"ijkl{ S»i (U^) - X(^) Sy-:r(u) + Sjj.y(ui) II (u-^) - x(^) Sfc;,^(u) + sw,j,(ui) 



/i(a;) < (x) dx 



X r 37 1 r X ' 

+ I 4jfci(2;)x(-) Sy,2;(u(x)) + Sy,j^(ui(x, -)) Sfei,^(u(x)) +Sfci,j^(ui(x, -)) 

n e L £ J L e . 

X r X 1 r X 

+ / dlj{x)x{-) di^x(p{x) + di^y(pi{x,-) dj^x(p{x) + dj^yifiix,-) 



dx 



X r X 

2 / C?-fe,(x)x(-) S^j (U^) SkLa:{u{x))+Skl,yiui{x,-)) 



dx 



dx 



2 / d'lJx)x{-) ip%x) d-j.xip{x) + dj^yipi{x,-) 



dx. 
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Using the ellipticity property of the elastic (cf^j,;) and dielectric (dfj) tensors, 
we get 

ac ||s,y (u^) - x(^)sy\x(u(a;)) - x{^)stj,y{ui{x, ^)) 
+ ad \\di ip^ix) - x{^)di^x'pix) - x{^)di^y(piix, ^) |j|2(o) 



< 



f.i{x) ul {x) dx 

4jkl{x)x{-)[Si3,x{^i^)) + '^^jA^li^^ -))J [Sfc;,x(u(x)) +SfcZ,y(Ui(x, -)) 

dij{x)xi-) di.x^ix) + di^yipi{x, -) dj,^ip{x) + dj^yifi{x, -) 
4jki{x)xi^) slj (u^) s ki, ^ iu{x)) + Ski,y iui{x,^)) 



dx 



dx 



dx 



dtjix)x{-) di ip^ix) dj^^ipix) + dj^yipi{x, -) 



dx. 



Using the fact that Sij^y{ui{x,y)) and di^yipi{x,y) are the admissible test func- 
tions and taking the limit in the sense of the two-scale convergence in the second 
right-hand side of the inequality, we obtain 



e->0 



Oic lim ||sy- (u^) - x{^){sij,x{u{x)) - Sij^y{ui{x, ^))} H^^^^^^ 



X, 



+ad lim (fi^ix) - x{-){d.i^xf{x) - d^^y^i{x, -)} 



< 



a Jv 



fi{x)ui{x) dx dy 



a Jy' 



Cijki{x,y)[s^j^^(u{x)) + Sij,y{ui{x,y))][skLx{u{x)) + Ski,y{ui{x,y))] dx dy 



dijix,y)[di^xp{x) + di^yipi{x,y)][dj,x(p{x) + dj,yipi{x,y)] dx dy (47) 

n Jv 

Recalling the form of the two-scale homogenized problem (fT0|) -(ITT |) - (fT2|) . we 
observe that the right-hand side of the inequality (|T7|) vanishes, so tha, we 
obtain 

lim \\sij (u^) - x( j){s»j,x(u(x)) - Sij,y{\ii{x, ^))} ||L2(a)= 



and 



X X 

lim \\di ip^{x) - xi-){di^x^ix) - di^y(pi{x, -)} \\mn)= 0. 
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6 Conclusion 



In this work, we have given the new convergence results, and the exphcite forms 
of the elastic, piezoelectric and dielectric homogenized coefficients. The two- 
scale convergence is applied to our problem yields the strong convergence result 
on the correctors. This technique of two-scale convergence can handle also other 
homogenization problems, in medium which has periodic structure for example 
the laminated piezocomposite materials or fiber materials (see [S] [5] [TT] [T3] 
|12p. Numerical implementation for perforated, laminated and fiber structures, 
will be presented in forthcoming publications (see [12j). 
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